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PREFACE 


This document Is submitted by the McDonnell Douglas Astronautics Company 
to the National Aeronautics and Space Administration and was prepared 
under Contract NASS- 331 91, "Statistical Energy Analysis of Complex Struc- 
tures." The study was directed by R. F. Davis. J. B. Herring of the 
Vibration Analysis Branch of the Systems Dynamics Laboratory of Marshall 
Space Flight Center administered and directed the contract. 
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SJWARY 


The aerodynamic and external acoustic environments of current aerospace ve- 
hicles generate very significant high-frequency random vibration structural 
response. This vibration response has proven to be a primary design con- 
sideration for the short life requirements (one flight) of past space vehicles. 
The effects of this adverse vibration environment will Increase for the 
reusable vehicles that are currently being developed for space exploration. 

The most efficient method to achieve optimum vehicle design for this high- 
frequency vibration environment Is to generate meaningful design and test 
criteria early In the design phase of the system and to periodically update 
these criteria throughout the various phases of vehicle development. The 
methods Included under the term "statistical energy analysis," or "SEA," 
provide a means of predicting such high-frequency vibration criteria In systems 
that do not conform well to other analysis methods. These energy analysis 
techniques are denoted "statistical" because they involve averaging structural 
responses over portions of the structure. This averaging is perl rmed over 
time and space and In frequency bands. 

Response predictions for a structure are made by modeling the structure as 
a number of elements, deriving power flow equations for each of these elements 
(Including acoustic or mechanical energy sources), and simultaneously solving 
the resultant system of equations for the element response levels. 

This report presents the results of an effort to document methods for accom- 
plishing such response predictions for craiwonly encountered aerospace 
structural configurations. The effort included application of these methods 
to specified aerospace structure to provide sample analyses. The report 
has been arranged in the form of an applications manual, with the structural 
analyses appended as example problems. Comparisons of the response 
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predictions with measured data are provided for three of the example 
problems. Other appendices provide a derivation of statistical energy 
analysis response equations* application guidelines* and response solution 
programs for programmable calculators. 
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SYMBOLS 


A 

a 

C/Cc 

eg 

Ct 

Co 

D 

E 

f 

6» 9 

h 

I 

L, I 

m 

N 

n(w) 

P 

r 

S 

t 

V 

V 

W, w 
w 

na 

^a,b 

K 

V 

ir 


area 

acceleration 

fraction of critical damping 

group velocity {- 1.07 [w Cj^ t]^) 

longitudinal wave velocity 

speed of sound In air 

dissipation of damping; bending stiffness 

total energy of an element, modulus of elasticity 

frequency 

gravitational acceleration 

thickness 

moment of inertia 

length 

mass 

number of modes 
modal density 
pressure 

radius of curvature 
power input 
thickness 
volume; velocity 
velocity 
width, weight 
weight density 
Incremental value 

damping loss factor = 2C/Cc] 

coupling loss factor (* ) 

radius of gyration 
Poisson's ratio 
3.14159 
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p 

a 

as 

Element 

System 

< > 


density 

radiation efficiency 

average mode-to-mode coupling between 
Elements a and b 

angular frequency, normally center frequency 
of a frequency band 


NOTATION 

a set of modes modeled as one unit of a system, 
all modes in a frequency band having identical 
energy (on the average) 

the total structure and associated energy 
sources under consideration (may be only a 
portion of an actual structure) 

indicates averaging over both time and space 


4 



INTRODUCTION 


Three methods are in coinnon usage for predicting the vibration response of 
aerospace structural systems: classical modal techniques, coitparative scaling 
using available data banks, and empirical formulas. Each of these nethods 
is limited in application by inherent characteristics of the method. Although 
classical dynamic analysis techniques for predicting dynamic response work 
well in the frequency range of the lower structural resonances, their appli- 
cation to high-frequency regimes is limited by model complexity, computer 
size capability, and cost, and is not amenable to rapid estimates. Use of 
data banks is limited to structural configurations which are very similar 
to the previous designs on which the data bank is based. Likewise, empirical 
formulas can, in general, be applied validly only to structures of the specific 
configuration for which they were derived and which match previous designs. 

The advent of reusable space vehicles featuring new configurations with unique 
forcing fields requires extending these techniques beyond the configurations 
from which they were developed. An alternative and relatively simple approach 
to nigh-frequency vibration analysis has been developed which is known as 
statistical energy analysis (SEA). 

Statistical energy methods have been developed to consider the distribution 
and transfer of energy among the modes of a vibrating system. These methods 
assume that the modes of a system being analyzed contain all the vibratory 
energy of that system. Therefore, for SEA to have valid application, all 
significant energy of a system must be "resonant" as opposed to "nonresonant." 
A parameter for evaluating this condition is examined in Reference 1. 

The SEA methods separate the frequency range of interest into frequency bands, 
which are analyzed independently. The methods assume that the energy in 
the modes of one frequency band is not transmitted (through coupling) to 
modes in other frequency bands, either within an element or among the ele- 
ments of a system. 
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An Important factor In validating the space and frequency averaging Inherent 
In SEA Is the number of modes Included In each frequency band. MIth many 
modes excited In one frequency band of an element « the vibratory energy may 
be expected to be well distributed throughout the element and among the 
various modes, and averaging will furnish a valid approximation of actual 
values. This effect is demonstrated In Figure 1 which shows, for a par 
cular structural syston, that predictions made with fewer than 20 modes p^r 
element per frequency band exhibit considerably more scatter than predictions 
with more contributing modes. Since constant percentage bandwidths (such 
as octave or one*th1rd*octave) are generally utilized to obto'i predictions, 
the narrower bandwidths In the lower frequencies result In fewer contributing 
inodes and therefore less accuracy at these lower frequencies. Accordingly, 

SEA Is generally applied to frequencies of 100 Hz or greater for typical 
aerospace structure. Also, model elements are chosen as generally gross 
portions of structure rather than representing fine details In order to 
maintain a high number of contrUuilng modes per element. An example of 
such modeling for a section of skirt structure on a launch vehicle would 
utilize three elements, one element representing the skIn/strInger external 
shell, another representing an equipment-mounting panel, and the third repre- 
senting the components on the panel (which are considered to be "smeared" 
over the panel In an average sense). 

The assumptions, then, upon which statistical energy analysis Is based are: 

A. The modes of the elements of a system contain all the '*^brator>' energy 
of the system. 

B. Only modes occurring within the same frequency band are coupled. 

C. The energy In one frequency band of a system element Is equally 
distributed among the modes of that element occurring In the frequency 
band. 

D. For two coupled elements, all of the modes occurring In one of the 
elements in one frequency band are equally coupled to each mode 
occurring in the same freqi zncy band In the other element. 
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DESCRIh^ ton 


The T'ollowtng section will demonstrate the application of SEA methods to a 
common structural systan of interest. While particular systans will require 
uh lue adaptations of the methods, the analytical procedure should be fairly 
ge..eral. 


Consider a section of airframe of an aerospace vehicle, such as a missile, 
reentry vehicle, or airplane, consisting of the external shell, an internally 
mounted equipment panel, and a ccnnponent mounted on the panel. 



Equipment 

Pml 




Figure 2. Typical Aerospace Equipment Panel Installation 

Several reasons may exist for analyzing the noted section of interest without 
performing an analysis that encoirpa^ses the entire vehicle: a local struc- 
tural change in the section for an operational vehicle, a design evaluation 
examining several locations and configurations for the panel, or evaluation 
of an alternate location for the component. 

This segment of structure can be represented with an SEA model of three 
elements as shown below. 


8 



Note that this iM)de1 will reduce to two elements for the ana ysis of a com> 
ponent mounted directly to the shell. 


Another ccmnon structural configuration of aerospace interest consists of a 
airframe section containing an internal bulkhead, with a coiqionent mounted 
on the bulkhead. 



Figure 3. Typical Aerospace Bulkhead Installation 

The corresponding SEA model for this structural system requires three elements 
and is of identical form to the shell/equipment panel/component model. 
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The analytical differences between the two strtictures will appear in the 
parameter values selected to represent the dynraiic properties of the various 
elements. There is no difference in the configuration of the two imxlels. 

The el&nents of these models are a very straightforward representatiwi of 
the structures. The only significant decision required is in definition of 
the amount of external skin to include with the model. Selection of the 
correct skin area yields a balanced system, with the energy flowing into the 
model subsystem mechanically fnun the renaint^r of the structure equal to the 
energy transported mechanically out of the subsystem. This situation therefore 
exhibits a zero net flow of energy across the boundaries of the model. In 
many cases this ideal condition can be approximately achieved by establishing 
the model boundaries at points halfway between major structural loading points, 
i.e., halfway between attach points of two adjacent equipment panels, halfway 
between the panel attach point and a fuel tank bulkhead, halfway between the 
panel attach point and a large component, etc. The effect on response pre- 
dictions of incorrect estimation of the skin area will be in essentially 
direct proportion to the error; selection of an area too large by 10% will 
result in predicted levels {g /Hz) that are too high by approximately 10%. 

Once the structural system has been modeled, the next step is to select the 
applicable equations for the model. These equations, listed below, were 
determined by examination (recognizing that thp coupling between the skin 
and conponent elanents is zero) of the SEA system equations listed in the 
Conclusions section of this report. 

(wHl + N 2(^.2 ) El - Ni 4)12E2 = Si 

-N 2 <j> 12 El +(un 2 Ni(j)i 2 + N3(|»23) E 2 ~ N2<|)23E3 = 0 

~N34>23E2 +(uri3 N24I23) E 3 = 0 


where 


0 ) = angular frequency (average) of system 

n, = Element e loss factor ( cVittcal ' Amplng ratio ) 

Ng * number of modes resonant in Element a 

b “ power transfer coefficient for coupling between modes through 
’ the structural joint (i|>« k * 'I'k ,) 

Cl f D D 
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Eg > total energy of Element a 

Sg > power Introduced Into Elenient a from an external source 

Each of these equations represents a power ^ow equation for one of ^e eloaents 
of the model. Together, the equations for. in algebraic system for the solution 
of the E{, provided the other terms can be evaluated for a structural system. 

The Input term. Si, will generally represent an acoustic excitation of the 
system. A derivation of the SEA response equations Is provided In Appendix I. 

The next step Is to evaluate modal density, daiiqpting, and coupling parameters 
for the system. 


Modal Density 

The shell elenent Is composed of the cylindrical skin, stringers and ring 
frames. The modal density of these subeleinents can be determined with the 
appropriate equations of Table 1 In Appendix II. 


Skll;: 


As / M 
AincpC^ Vwr-/ 


Stringers (assimilng plate-type response for high frequencies): 




Ring frames (assuming beam- type response): 




I 1 _ 


where 

N 

Og = modal density of Element a “ ^ 

iuu - frequency bandwidth selected for analysis 

As = surface area 

icp = radius of gyration of plate cross section 
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Cl • longitudinal wave velocity 

c* 

ay * structure ring frequency » 
r » radius of curvature 
L * length 

icb ~ radius of gyration of beam 

Sunmlng these contributions. 

"> • "skin * "str * "rf 
and 

Nj » nj(^) 


The equipment panel (or alternately, the bulkhead) Is a ribbed plate which 
can be subdividrd Into plates and beams: 

bulkhead “ *'plate * *^r1bs 

The component can generally also be subdivided Into plate-, beam-, or shaft- 
type elements, depending upon specific design. If, however, experimental mc^al 
response data are available on similar compo.-ents, an estimate of the com- 
ponent modal density can be obtained from a plot of mode number versus f. equency. 
The slope of the approximate line joining the points, as indicated in fig- 
ure 4, Is the modal density. 
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Da ncing 

Th :2 dairping paraneter values for the elements should be estimated based on 
experience with similar structures. References are provided In Appendix II 
which win isslst the user In selection of values. 

Couc l ng 

Appenoix I. provides coupling factors that apply for a number of structural 
joints, tor the typical structure under consideration, coupling factors are 
p-o' ded for specific skin/equipment panel joint configurations (also for 
SNin/bulkhead joint configurations of the alternate structural system) that 
are similar to many aerospace installations. 

A wide Variety of configurations may be encountered for the component/equ1[»nent 
panel joint. The first approach to evaluating coupling factors for this 
Joint is to check Appendix II and other available sources for a similar joint. 
When coupling factors for a similar joint cannot be found, this parameter 
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may be determined through an SEA evaluation of response data with a similar 
joint. The SEA parameter for the similar joint may be evaluated by the 
following procedure. 

The response relationship between a component and mounting panel 1s» In 
general, p.'ovided by an equation In the form of the third of the set of 
SEA equations: 

"Nc*!>pc^p '*■ Np(j>pc) Eg = 0 

where the subscript notation is: c = conq)onent 

p = equipment panel 

This equation provides a means of solving for the coupling factor If the 
other system parameters can be evaluated and response data are available. 




The total energy terra for an element Is provided by 
~r <a^i> 

El * roi <Vj*> = <"1 -jjz- 


m, = element i mass 

V,* = element i velocity 

aj = element i acceleration 

< > indicates averaging over time and over area 

which yields 
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This equation defines coupling factor based on the relative response 
level, , of the equipment panel to the component for the similar joint 

<a*g> 

configuration. (Similarity to the prlnairy structural system under analysis 
assists In defining Che SEA parameters In the equation, generally reducing 
the effort required to determine the coupling parameter using resptmse data.) 

The final step before obtaining the response solution for the system Is to 
define Si, the Input term. For the sample structure, excited by a reverberant 
acoustic field, this term Is 

Sj * ^ (surface) 

ojo v'aw) ra-f 


Co » speed of sound In fluid medium 
» surface area of Elen«nt 1 
P * acoustic pressure 
a » radiation efficiency 

s number of surface modes of Element i (excludes modes of ring 
frames, stiffeners, etc.) 

m<j =« mass of Element i 

Appendix II provides radiation efficiency values for both flat panels and 
circular cylinders. 

The parameter values ma>' now be substituted into the system of SEA equations 
(for each bandwidth of interest) and the response solution for the Ci obtained. 
The previously noted relation. 



can then be used to present the element responses In the form of 




(for response in g's) or 7-- . --r (in g^/Hz), where g is the gravitational 

(Auig*" 


acceleration constant. 
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SUMMARY OF STATISTICAL ENERGY ANALYSIS 
LIMITATIONS AND EQUATIONS 


The nethods presented In this document will provide estimates of the high- 
frequency vibration environment for structural systems. The user should 
be aware of the requirements and limitations for the application of these 
methods. A summary of the principal limitations Is provided below. 

1. Application Is more valid in frequency ranges where many modes are 
excited. Care must be taken In evaluating the lower frequency limit 
of applicability for a particular structural system. 

2. Application Is valid only for systems containing all their energy In 
modal resonances, therefore SEA ^''^s not apply to heavily damped systems. 
Reference 1 provides a means of evaluating this requirement for specific 
systems. 

3. Response predictions determined with these methods represent averages 
over generally gross portions of structure. Therefore caution should be 
taken In applying these average values with nonuni formly configured 
structural elements such as panels with relatively massive Integral 
stiffeners, where response amplitude of the panel segments mey be expected 
to differ considerably from that on the stiffeners. 

A summary of the SEA response prediction equations is provided below. In- 
spection of the system equations will indicate the terms required to expand 
the set of equations to acconmodate a system with any number of elements. 
Likewise, the simplification possible for systens which do not have each 
element connected to every one of the other elements can be determined by 
setting the appropriate (j)jj terms to zero. 

Guidelines for structural modeling and parameter evaluation are provided 
in Appendix II. 
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Statistical Enemy Analysis Equations for a Four-Element Systaw 


(wHl + N2V12 Ns(j)l 3 + N 4 lj>l 4 ) El - Ni4>12E2 - Ni<|> 13 E 3 - Ni(i>i 4 E 4 * Si 

“N2<|>12Ei +(ujri2 + Ni(|>i 2 N3<j>23 N4<|)24 ) E2 “ N24>23Es ~ N2<^2*»E4 = S2 

“N 3 ^ 13 Ei - N3^2sE 2 ■•■(wns "*■ Ni(j>i3 + N 24>23 N4({)34 ) E3 - N34>34E4 * $3 

“N44I14E1 ” N 44 > 24 E 2 “ N4(|)34E3 +{wn 4 Ni<j)i4 + N 2^24 ^3(1)34 ) E4 ® S4 


* angular frequency (center frequency of analysis band) 

Hj = element damping factor 

N{ = number of modes excited In the element In frequency band of analysis 
<l*ij = element coupling value (symnetric: (j)jj = 



mi 


<a"i> 




angular acceleration 

. ^ ^ Zir^Cn^Ai <Pi^>oi . . .. 

energy Input term = -a J- Ni (for reverberant acoustic 

o)q (Aw) mj excitation only) 

local speed of sound In surrounding medium 
acoustically excited surface area of element 
acoustic pressure 
radiation efficiency 
mass of element 


These equations are frequently expressed In an alternate format by combining 
terms Into a coupling loss factor: 


„ - 
'>ij - » 


= coupling loss factor (nonsymmetric; 


nij nji) 
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Appendix I 

DERIVATION OF STATISTICAL ENERGY ANALYSIS RESPONSE EQUATIONS 


Consider a simple structure modeled as the two-element system in the following 
schematic. 



In this schematic of the two elanents, a and b, the following nomenclature 
Is used: 

Sg = power Introduced Into Element a from an external source. 

Og = power dissipated within Element a. 

Pa^b “ power transmitted from Element a to Element b (= -Pb,a)* 

These values and the following derivations are for only a single frequency 
band; solution for the complete spectrum of Interest is accomplished by 
summing the predictions for the contributing frequency bands. 

Power flow equations for all of the energy passing through the two elements 
may b, expressed as 


Db - Pa.b * h 
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The energy dissipated per unit time Is defined in terras of the element loss 
factor as 

h = “HaEa 

where 

(0 - angular frequency (average) of system 
% = Element a lots factor ( critlc T rdLnp T n ? ratfo ~) 

Ea = total energy of Element a 

The net power transmitted from the resonant modes of Elanent a resoiant 

modes of Element b is 


Pa,b = Va,b^a - Na(|)a,bEb 

* (power transmitted from b to a) - (power transmitted from a to b) 


where 

Na - number of modes resonant in Element a 


ij>a b ® power transfer coefficient for coupling between modes through 
the structural joint (<|)a,b “ ^b,a^ 


Performing the Indicated substitutions, the power flow equations become 

wHaEa Nb4^a,bEa " Na<|)a,bEb “ 
tonbEb + Na1>a,bEb " '^b‘l’a,b^a ^ ^b 

For systems with the parameters n, N, (j> defined and the inputs, S, known, 
the power flow equations form a set of linear, simultaneous equations for 
the unknowns Ea and Ep at the average frequency w. 
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App&idix II 

GUIDELINES FOR APPLICATION OF STATISTI** ENERGY ANALYSIS 


PeflRltion of Structural Models 

One of the Initial steps In the SEA applications procedure will be selection 
and definition of suitable models. The basic considerations o^ nrdaling 
are to (1) determine the structural definition and detail requir.J, (2) evalu- 
ate energy sources, and (3) partition the significant portion of the structure 
Into the actual model elements. The first consideration, requirements for 
structural d«>f1n1t1on. Is to Insure that the model will both provide the 
Infonr'tlon desired and omit useless details. Structural assenblles may 
be lunq)ed together as a single element If finer definition Is not required 
since SEA uses averaged quantities, and averaging Is equally valid for mul- 
tiple portions of a structure as for a single part. Such lunq)1ng of elements 
also reduces tne bookkeeping associated with the analysis. 

The second consideration In SEA modeling Ir ev^^luatlon of the energy sources. 
This consideration assists In limiting the size of a model. Basically, any 
structural boundary across which the net energy flow Is zero reoresents a 
limit to the need for modeling. 

The final step In modeling Is to actually partition the significant structure 
Into elements In line .th the previously stated principles. The elements 
* represent generally gross, continuous portions of the structure. 

Damping 

The structural damping must be defined for each element of the models as 
one of the procedural steps. This parameter Is not unique to SEA and must 
appear In some form In every response analysis. While much Investigation 
of structural damping has been ccompllshed, and over a l^ng period of time, 
selection of appropriate values remains very much a matter of engineering 
judgement based on past experience. References 1 and 2 provide Information 
useful In defining the damping of st."uctures. 
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Modal Density 


Modal density Is the para>neter which is used to evaluate the number of reso- 
nant inodes present within a particular frequency band of a given structural 
subset. Approximation equations are presented in Table Il-l (from Inference 4) 
rdiich define this parameter for specific structural shapes. 

Alternate nethods of evaluating the modal densities of specific structural 
configurations are available. One alternate method makes use of computer 
programs generally available for analyzing the response of commonly encoun- 
tered structural shapes such as pinned-end cylinders, liquid-filled cylinders, 
etc. These pn^rams are utilized to determine the frequency respwise of 
the elements, thus yielding directly the number of resonant modes within 
the frequency band. 

A second alternate method has utilized classical low-frequency modal analyses 
of the structure to evaluate the modal density for model elements. This 
nethod, of course, is only practical for structures which have previously 
received modal analyses, or if a low-frequency modal analysis is being per- 
formed in conjunction with the high-frequency SEA prediction. This method, 
indicated in Figure II-l, involves graphically plotting the response frequen- 
cies from the modal analysis versus mode number for the structure which is 
being represented by an SEA model element. The slope of the plotted points 
can then be determined which, assuming extension of the curve into the high 
frequencies to be valid, yields the value for modal density of the element. 

Structural Coupling 

The SEA parameter for the structural coupling between elements is unique 
to this form of analysis, and its definition represents one of the most 
significant steps in the procedure. As a consequence of the relative 
newness and lack of previous application of the SEA methods, very little 
in^^ormation has been available concerning appi^opridte vali es of this parame- 
ter for structural joints in general. The coupling values for tv/o typical 
joints are shown in Figures II-2 and n ' 



Table II-l 

MODAL DENSITIES OF SOME UNIFORM SYSTEMS 
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■^See next for definitions of symbols. 


Symbol Definitions for Table 1 1-1 



cross- section area 
surface area 

acoustic wave velocity 

Longitudinal wave velocity 

membrane wave velocity 

string wave velocity 

torsional wave velocity 

plate rigidity 
Young's modulus 
shear modulus 
thickness 

centroidal moment of inertia of A 
polar moment of inertia of A 
torsional constant of A 
length 

membrane tension force/unit edge length 

string tension force 

volume 

radius of gyration of A 

radius of gyration of plate cross section 

Poisson's ratio 
frequency (radians/time) 
material density 


OOt/€U.A|i 

o 
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Figure II-l. Approximation of Nodal Density by Graphical Nethid (Delta 
Fairing Nodes) 




Figure II-2. SEA Coupling Parameter for Small Reentry and Intercept 

Vehicle Field Joints (Typical, from UpSTAGE Ground Test Data) 
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Figure II-3. SEA Coupling Paraoneter for Launch Vehicle Tank/Skirt Joint 
(T^ical, Based on S-II Data) 


Reference 6 provides equations for evaluating the coupling values of various 
beam and plate joints. Two of the nost useful of these relations are pre- 
sented below. 

Two plates of approximately equal stiffness joined at right angles: 


<Pl2 


c 


9 


ttA, Hz 27 
1.07(o)Cj t)^* 


L = joint length 
t = plate thickness 
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Beam cantilevered to a plate of equal thickness: 


. . 2irf W 

Np - nimber of inodes in plate 
M = width of beam 
t = length of beam 


Acoustic Coupling 

The input term will generally involve a transfer function to couple a fluc- 
tuating pressure field to the structural system. A reverberant acoustic 
field may be coupled to a structure with the relation presented in the Con- 
clusions section of this report. Use of this expression for predicting 
response from other acoustic fields requires the definition of an "equivalent" 
reverberant field, or the coupling terms must be modified. The developirent 
of the reverberant coupling terms will Indicate an approach that could be 
used in defining coupling terms for other pressure fields. 

The radiation efficiency term, o, appearing in the input relation for re- 
verberant acoustic fields of the Conclusions section may be determined 
from Figures II-4 and II-5, which are taken from Reference 5. 
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TMfRO OCTAVC 0AMO CE.^TC« FRCOUCNCt IN CYCCtS PC« SfCOWO 


Figure II-5. Radiation Efficiency o of a Cylindrical Shell (The peak in 
the radiation efficiency about the ring frequency f^ is as- 
sociated with increase in wavei.peed due to curvature.) 

(From Reference 6) 
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Appendix III 
EXAMPLE PROBLEMS 
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Example Problem Number 1 

SPACE SHUHLE EXTERNAL TANK - UNLOADED STRUCTURE 

The structure to be analyzed is located in the Space Shuttle Extr.nal Tank 
intertank area at 270® (-Y) on the station 1034.2 frame. This location 
corresponds to a vibration measurement location in use during Main Propul- 
sion Test Article (MPTA) testing for the Space Shuttle. The measurement 
location is indicated in Figure III-l. This location is on a ring frame 
which is surrounded by skin panels with external stringers. This portion 
of the structure is not loaded by component installations and can therefore 
be considered typical of unloaded aerospace shell structures. 

Model 

The unloaded structure can be represented by the simplest of SEA models, 
consisting of a si. le element excited by the external acoustic field. 

The model element will include the area 45° to either side of the measure- 
ment location (225° to 315°) and halfway to the adjacent frames at stations 
985 and 1082. This element is indicated in Figures III-2 and III-3. 

Response Equation 

The SEA response equation for the one-element system is 
wii El = Si 

This equation simply equates the energy dissipated by damping within the 
element to the energy transmitted from the external acoustic field. 
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Figure III-l. MPTA Vibration Measurement Locations (Station 1034) 
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Figure I I 1-3. SEA Model Element for Unloaded Structure 
(Cross Section) 
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Damping 

The damping parameter was evaluated with response data measured during an 
acoustic fatigue test of Saturn S-IV8/V Interstage panels (Reference III-l). 
These panels feature a skIn/stringer construction similar to the current 
structure and are cf similar size. The approximate relation for damping 
of n = where (Af)® represents the bandwidth to the half-power points, 

was evaluated for a number of response measurements on the test. The 
resulting values have been plotted in Figure III-4. Straight lines have 
been faired through these data for a simple graphical approximation to 
the damping which Is used for this analysis. The approximation lines 
have been positioned on the low side of the obvious mean of the data since 
low values for damping result In high predicted responses that are con- 
servative for design purposes. 

Element Energy 

The energy of the model element Is represented by 

Ex = mi<Vx > » mx — ^ 

( 0 ^ 

The elanent mass was estimated from available detail drawings by summing 
the volume of the element subsections (skin panels, stringers, frame sec- 
tions, etc.) and multiplying by the material density. 

310 lbs 

g— 

Acoustic Power Input 

The external acoustic field Is assumed to be reverberant, therefore the 
input term can be represented by 

. (surface) 

1 <- 

Wq" (Ao)) mj 
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selected for current analysis) 



The surface of the subject structure Is composed partly of skin panels 
and partly of stringers. The correct representation of the term 
Is thus ' 

^jp ^ ^:s ^IS 
n>ip n»is 

where the subscripts p and s Indicate panel and stringer values, respec- 
tively. Since the stringers tend to bound the surface Into plate areas 
of relatively small curvature, the appraxlmate relation for high frequency 
modal density of plates, 

"(“) ■ 

was utilized for both skin and stringer arebs. 

f 

N « n(Aw) 

also m = ^ V * At 

w = weight density (7075 A1 * .101 Ib/in*) 

V = volume 
t » thic^ness 


yielding 
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The surface area of the element is approximately 48.5 in x 259 In, with 
36 stringers of the cross-sectional dimensions indicated below. 





The effective width of the stringer is approximately 8.5 inches. 

As = 36(8.5 X 48.5) = 14,841 in^ 

The skin panels are 0.071 inch thick and have an area of 
Ip = 48.5 X [258 - 36(4.41)] = 4862 in^ 

Then 

AN - ^f Jl 2g(U^ ' . ^s] 

2 1 V7E LV VJ 

= 151,934 Af 


Acoustic levels me. sured during the MPTA test are presented in Figure III-5. 
The three measurements were averaged for each one-third octave band center 
frequency and the average value used as the required acoustic pressure input 


<P^> = 10 exp 


(SPL) 


avq 


10 


X 8.41 X 10 


■10 


The values for <P^-> are listed in Table riI-1. 
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SOUitxJ pressure level - tiB (re 2 X 10*'' H/M 



0NE-TBI8D OCTAVE BMID FREQUENCIES > HESIE 


Figure II 1-5. Space Shuttle External Tank Main Propulsion Test Acoustic tota. 

Three measurements external to intertank area (solid lines); 
upper, dashed line represents measurement subsequent to analysis 
(see text). 
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Table III-l 


f 

Ih. 

<P*> 

IM 

o 

50 

0.52 

116.5 

.000794 

63 

0.&3 

117.2 

.00100 

80 

0.48 

114.4 

.00126 

100 

0.47 

116.6 

.00158 

125 

0.46 

116.6 

.00200 

160 

0.45 

118.1 

.00251 

200 

0.44 

118.9 

.00316 

250 

0.3b 

118.9 

.00398 

315 

0.25 

120.0 

.00631 

400 

0.20 

122.2 

.00316 

500 

0.15 

121.7 

.00251 

630 

0.12 

121.9 

.00251 

800 

0.090 

120.8 

.00251 

1000 

0.070 

118.3 

.0100 

1250 

0.054 

116.8 

.0200 

1600 

0.040 

115.7 

.0316 

20CO 

0.032 

114.7 

.0631 
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Radiation efficiency values are taken from Appendix II, based on a coinci- 
dence frequency of = 7683 Hz (for the 0.063-inch thickness of the 
stringers since they furnish most of the surface area) and a ring frequency 
of fp = 198 Hz. These values are also listed in Table III-l. 

Response Solution 

The response predictions for the element were determined in each one-third 
octave band from 50 to 2000 Hz. A sample prediction for tte 50 Hz octave 
band Is presented below as an example. 


wHi ~ Sj 


^n"[(1116 X 12)^] [(151 ,934) 8.41 x lo'**] [■0007?x 3 

[(2x X Mj-'jfih, X 


where Af = for one-third octave bands. 

Dividing by g and solving for the mean squared acceleration: 
<a.*> 

- 4 - = 0.0130 

r 

The acceleration spectral density level is 

<ai*> _ .0130 

g“Mf) ' /J0.\ ' 

I 4.33/ 


and the root-mean-squared acceleration in this one-third octave band is 



0.114 
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The predicted response levels are plotted in Figures III-6 and III-7 for 
the respective acceleration spectral density and 9^03 in one-third octave 
bands. 

The increasing levels above 1000 Hz which are most noticeable in Figure III-7 
are an unexpected result. Inspection of the input parameter values for 
the response solution shows this result can be attributed to the radiation 
efficiency. The value for this parameter is controlled by the coincidence 
frequency of f^ = 7683 Hz for the 0.063- inch-thick stringers. This fre- 
quency corresponds to the peak radiation efficiency values and is higher 
than typically encountered with aerospace structures, thus causing the 
radiation efficiency to be still increasing at 2000 Hz with the resulting 
high predicted levels. 

Reference I I 1-4 danonstrates that ±3 dB accuracy may be expected for an 
SEA response prediction when 20 or modes per analysis band are excited 
in each element or above the ring frequency for cylindrical structure. 

That result was based on a relatively small, s^iff vehicle with an ellip- 
tical shape. These requirements correspond to 20 Hz (20 modes) and 200 Hz 
(ring frequency), respectively, for the current structure. Therefore, 
the predictions may be expected to demonstrate ±3 dB accuracy at frequencies 
above 200 Hz, and probably above 20 Hz. 

Comparison of Prediction wi -h Measured Test Data 

Subsequent to completion of this prediction, MPTA response data were made 
available for comparison. However, one of the acoustic measurements pro- 
vided with the response data showed an increase of more than 10 dB throughout 
the spectrum (see Figure III-5), while another measurement agreed with the 
previous acoustic data. Further investigation revealed that the initial 
acoustic data were for a 70% thrust level, and that the high measured 
levels were probably valid for one side (adjacent to Orbiter engines) 
of the external tan!, at the 100% thrust level with a reduction in acoustic 
levels around the tank to the opposite side. Because of the resulting 
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Acceleration Spectral Density, g^/Hz 



Figure I1I-6. Measured Test Data (dashed line) vs. Predicted Response for 
Two Acoustic Input Levels (solid lines) on External Tank 
Interstage Area - Unloaded Shell Structure 
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uncertainty in the exact input acoustic levels, the response predictions 
for the data comparison are presented in Figure III-6 for two input levels 
the original, averaged value and the high measurement level - with the 
expectation that the correct level actually lies oetween the two. The 
vibration response data correspond with this expectation and lie chiefly 
between the two predictions (SEA response predictions should be compared 
to average values of the response over one- third octave bands rather than 
peak response values). At 1000 Hz, the predicted response exhibits a 
change in slope and begins to increase at higher frequencies, a result 
determined to be due to an increase in radiation efficiency values near 
the coincidence frequency of the external stringers. This discrepancy 
in response with the test data is most likely due to improper definition 
of damping values about the coincidence frequency. This overprediction 
of response relative to the measured levels would lead to a conservative 
result when used for design purposes. 
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Example Problem Number 2 
SPACE SHUTTLE EXTERNAL TANK - LOADED STRUCTURE 


The structure to be analyzed Is located In the Space Shuttle External Tank 
Intertank area at about 200 ° on the station 1034.2 frame (refer to Figure III>1). 
This location corresponds to a vibration measurement location used during 
Main Propulsion Test Article (MPTA) testing for the Space Shuttle, as for 
Example Problem 1. The structure is identical to the skin/stringer/ring 
frame structure of Example Problem Number 1. but is loaded by the 260 pound 
DFI box installation and can therefore be considered as typical for aerospace 
shell structure loaded by heavy conqionents. 

Model 

The configuration to be analyzed represents structure which is loaded by the 
OFI package. This package mounts on supports between frames. The model for 
this structure will have two elements, consisting of the external shell and 
the OFI package with support intercostal s. The shell structure to be included 
will extend halfway to the fran«s adjacent to those carrying the support 
structure and a circumferential width identical to twice the support frame 
width, as Indicated by Figures III-8 and III-9. 

Response Equations 

The SEA response equations for the two-element system with external acoustic 
excitation are 


(wTll + N2<}»12 ) El - Nl'^l^E2 * Si 

-N 2 <^ 12 Ei + ((jiT\z + Ni<J)i 2 ) £2 = 0 


where the subscript 1 denotes external shell values, and 2 denotes DFI box 
and intercostal values. 
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station 



Figure III-8. Structural Configuration for SEA Model of Loaded Structure 


Station 


Area of SEA Element Support Brackets 



Figure I I 1-9. Extent of SEA Model for Loaded Structure 
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Dffloplnq 

The damping parameter for the external shell Is Identical to that of Example 
Problem 1 in Figure III-4. 

Test experience with smaller electronic packages during the Delta program 
Indicates respo.ise amplifications of Q» 6 to 10 should be expected for the 
DFI package. Therefore, a value for the loss factor of n = 0.1 (Q*10) was 
adopted for this analysis. 

Modal Density 

The expression for number of modes, N{, required In the response equation Is 
determined by 

Ni = n| (Af) 


where 

n^ » modal density 

Af = bandwidth of analysis 


The portion of Intertank structure which has been designated as the external 
shell element Is actually composed of several hundred Individual parts (skin 
panels, stringers, stiffeners, ring frame segments, fittings, brackets, etc.). 
The modeling of all these parts by one element Is a feature of the averaging 
assumptions of the SEA approach and Is one of the most attractive aspects of 
SEA. Almost without exception, the Individual parts are plates or are formed 
with multiple plate sections. Therefore, the shell element modal density was 
calculated by simmlng the modal density of the Individual plate sections de- 
termined by the approximate relation for high frequency modal density of plates. 


ni(f) = I- 


2 Kp 


1 


2 Eg 

V12W(l-v*) 



» 6.52 modes/Hz 
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A 

where I ^ represents the sumnation of val'jes for all parts and their subsec- 
tions. 


The second element is made up of the DPI package plus the supporting inter- 
costal structure. The DPI package consists of a box, whose plate element 
modal densities can be determined as above, plus a panel loaded with electronic 
coioponents. Reference 2 indicates that the loaded panel will exhibit a greater 
stiffness (arJ resultant lower modal density) than an identical unloaded panel. 
An increase in stiffness by a factor of 2 was assumed for the loaded panel. 

The resulting modal density For the boj( element is 

noFiCf) = l"2)box * ^^^Vanel 

A] + _L f A\ 

^'box /T ^ panels 

\h2w(l 

* 0.27 + 0,12 
= 0.39 modes/Hz 



The modal density for the supporting intercostals s also determined with 
the plate equation: 


nj(f) = 




Eg 


(l-v^) 



0.32 


Therefore the total modal density for this element Is 
na(f) = nppj(f) + nj{f) = 0.71 

Modal Coupling 

The model elements are coupled by the joint between the vehicle shell and the 
DPI support intercostals. This joint is essentially two plates joined at 
right angles. Por the case of equal plate stiffnesses which is approximately 
satisfied (.071 inch Intercostals and .071 inch skin with supplementary stiff- 
eners), Reference 3 gives the relation for coupling loss factor of 
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TIj2 


Cg L I±] 

2ir*fAi \27/ 


Cg * 1 •07(b) Cj^ t) 


Vi 


orr 

J wO -V*) 


C)6 

L * Joint length 


From the basic def v ition 


hl2 


<}»12N2 


b) 




SO that 


♦12 


, _bL -g-qL / 8 \ ^ IJO 
N 2 2iT^fAi \27/ 


which makes use 0 ^ “ 4.33 for the one-third octave bandwidths to be used 

for analysis. 


Element Energy 

The element energy was handled as In Example Problem 1, 


c - m ^ 

^1 “ "'I ~zr 


b)* 


mi = 


484 lbs 


m 2 = 


^ 280 lb s 
9 
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Acoustic Power Input 


This teim was also handled as in Example Problem 1. 



The applicable term for 


(Mi) 

' 'surface 


237.473 Af 


R esponse Solution 

The response predictions for each element were determined in each one-third 
octave band from BO to 2000 Hz. A sample prediction for the 50 Hz band is 
presented below as an example: 

(wrii Nacbia ) Ei - Ni<j>i2t2 * Si 
-N2(|)12Ei + (wTl2 + Ni4>12) E 2 ® 0 


Substituting ^'or the parameter values, the expressions become 


X 50]t0.5J] . 


■ [®-®K4.33)][i?] 

2ir2[(1116 X 121*! 

290 

< ^2 ^ 

][l0^ 8.41 X 10"**] 

[.000794] 

.9 (2ir X 50)J 

1 [(233.473) 


(2tr X 50)*] 

^ 4.33] 
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and 


♦ jc2. X 50]C0.n . [6.52{^)][^]j[f - 0 


where Af = —r for one-third octave bands. 
4.33 


Dividing by g and solving for the mean squared accelerations: 

<a,*> dii 

“^2— =0.0130 -jr = 0.000662 


The acceleration spectral density levels are 


^ 0.0130 
9^^ / 50 \ 

U.33/ 


0.00113, 



g^Af) 


0.000662 



0.0000572 


and the root-maan-squared acceleration in this one-third octave band is 


(9™s>. 0.114 

» ^.000662" = 0.0257 

The predicted response levels are plotted in Figures III-IO through 111-12 
in acceleration spectral density and in one-t.iird octave band formats. 

The criterion of 20 modes per analysis band iti each elanent as a requirement 
for SEA prediction accuracy is satisfied for the moael at 125 Hz, indicating 
predictions with t3 dB occuracy above this frequency. 
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Acceleration Spectral Density, g^/Hz 



FREQUENCY 4HZ) 


Figure III-IO. Predicted Response (g*/Hz) for External Tank Interstage Area 
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Figure III-ll. Predicted Response (gnns) External Tank Interstage Area - Loaded Shell Structure 

(Element i) 






BAND NUMSenS 
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Figure 1 1 1- 12. Predicted Response (grms) External Tank Interstage Area - DFI Box and Supports 
(Element 2) 





Comparison of Prediction with Measured Test Data 

The structures analyzed in Example Problems Numbers 1 and 2 were selected 
to provide an evaluation of SEA methods in application to loaded and un- 
loaded structure. To circumvent i . effect of the uncertainty in the 
acoustic input levels as noted in Example Problan Nunber I, the available 
test data have been compared to the predicted values in a relative sense 
between the loaded structure of this example problem and the unloaded 
structure of Example Problem Nunber 1. A comparison of the element 1 
response in Figure III-IO with the lower predicted curve of Figure III-6 
shove the structural loading to have no effect on the predicted values 
which are essentially identical. Figure III-13 shows the approximate 
averii'T spect«*un values {faired graphically through the data) for the 
test measurements. These data show the loaded structure to have reduced 
response relative to the unloaded structure below approximately 800 Hz, 
and essentially the same average spectrum levels above that frequency. 

For this case and based on these specific measurement locations, the SEA 
method furnishes accuracy within 3 dB only above 500 Hz, which corresponds 
to 80 modes per analysis band in each element. 


Prediction of the relative response of the loaded structure through tradi- 
tional mass scaling would result in a reduction of response by the ratio 


434 lbs 

484 lbs + 280 lbs 


0.630 


or approximately 2 dB. The SEA method obviously yields more accurate 
results in the higher frequencies. 
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Relative Acceleration Spectral Density* dB 



10 100 1000 I'.OOO 
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Figure III-13. Comparison of Relative Response for Test Measurements on 

External Tank Interstage Area Unloaded and Loaded Structure - 
Approximate Average Values 


."jO U>6 1 16 JU^C 711 
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Example Problen Number 3 

SPACE SHUTTLE EXTERNAL TANK ~ DETAILED ANALYSIS OF LOADED STRUCTURE 

The structure to be analyzed is identical to that of Example Problem Number 2 
and is located in the Space Shuttle External Tank intertank area at about 
200“ on the station 1034.2 frame (refer to Figure III-l). This location 
corresponds to vibration measurement locations used during Main Propulsion 
Test Article (MPTA) testing for the Space Shuttle, as for Example Problems 1 
and 2. The structure is identical to the skin/stringer/ring frame structure 
with DFI box installation of Example Problem Number 2, but will be modeled 
in more detail, including internal acoustic excitation, to provide a response 
comparison for a measurement location on an equipment oanel inside the DFI box. 

Model 

The configuration to be analyzed represents structure which is loaded by the 
DFI package. This package mounts on supports between frames. The model for 
this structure will have four elements, corio'~t1ng of external shell, DFI 
support intercostals, DFI box equipment panel, and DFI box cover. The shell 
structure to be included will extend halfway to the frames adjacent to those 
carrying the support structure and a circumferential width identical to 
twice the support frame width, as indicated by Figures III-14 and III-15. 

Response Equations 

The SEA response equations for the four-element system with external acoustic 
excitation are: 


(ojrii N24 >i 2) El - Ni4’i 2E2 = Si 
-N 24 ' 12 Ei ■*'(»jjn2 Ni<t>i2 + N 3 <{> 23 ) E2 ~ N 2 <t 23 E 3 = 0 

~N 34 > 23 E 2 ■*"(wr)3 N 2 'J *23 Nu^^ai,) E3 = N34*34 E4 = 0 

-N(,4*3i*E 3 + (wn** N3<*>34 ) E4 = S4 
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Figure III-14. Structural Configuration for SEA Model 


Area of SEA Elanent 


Support Brackets 


Station 



Figure III-15. Extent of SEA Model for Loaded Structure 
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where subscript 1 denotes external shell values, 2 denotes intercostal values, 
3 denotes OFI box equipment panel values, and 4 denotes DFI box cover values. 

Damping 

The damping parameter for the external shell is identical to that of Exonple 
Problem 1 in Figure III-4. 

The dandling parameter for the DFI box/intercostal elanent of Example Problem 2 
was assigned a loss factor value of n = 0.1. This value was also adopted 
in the current example problem for each of the three elements (panel, cover, 
intercostals) resulting from subdivision of the previous DFI box/intercostal 
element. 

Modal Density 

The modal density parameters were calculated during Example Problem 2 as 

ni(f) = 6.52 modes/Hz 
nz(f) = 0.32 modes/Hz 
03(f) = 0.12 modes/Hz 
04(f) = 0.27 modes/Hz 


Modal Coupling 

Three joints are inck .eu with the model. The shell/intercostal joint was 
evaluatec* during Example Problem 2 and assigned the value 



V I 


The DFI box i;. connected to the intercostals by eight flanges which connect 
to right-angle intercostal brackets. Therefore the joints were considered 
as two plates jointed at right angles and evaluated following the approach 
used for the shell /intercostal joint. 

'^9'- (A.) = 2M 

Ni A, V27/ /f 
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The joint betti^en the DPI box cover and equipment panel Is also a right 
angle connection between two plates and was evaluated as In Exanqile Prob* 
lem 2. 


(t) 


CgL /JL\ , 255 
K 2iT*fA3 V27/ /P 


The effective thickness of the laminated aluminum and balsa plates required 
to evaluate the term 

Cg » 1.07(u»C£t)*''^ 


was determined through evaluation of the bending rigidity, D « mi4 Cj , 
for an equivalent single-layer plate of aluminum 

|< Lj >1 

Aluminuro/Balsa Plate: hi ( h 2 Balsa 

Jl-L*::. - 


M 


The equivalent laminated all-aluminum plate with the same rigidity has the 
dimensions |<- 

I ^ 

hi h2 1 ^ — L 2 

i_.ci= 


1 ; 


where l 2 * ■ 0.04 Li, since In the center layer 

2 ^ 
Dbalsa “ "»B*''cJg - ^PbVb)t|(^)= ^ h ‘ * U • 1) 


where the volume term, V, is evaluated per unit depth. 
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For the equivalent center section of aluminum: 


12 ^ At * (^2 •‘-2 • 1 ) 13 ^At 


Equating the terms 


•^balsa “ *^At 


ha 


(ha * Li • 1 ) ^ Eg “ (f *2 • *-2 * 1 ^ ^At 


or 




The effective thickness of a single-layer plate of alioninum vrith the same 
rigidity is defined by 

% “ ^^At 

or 

nUp 

<2 = ( k ')2 


For the equivalent laminated aluminum plate. 



and for a single-layer plate. 


yielding t^ * h? - .96 h| 
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For the DFI box equipment panel, hj * 0.96 in, ha » 0.92 in, so 
t = 0.33 in 

is the required effective thickness. 


Element Energy 

The element energy was handled as in Example Problems 1 and 2: 

_ _ 4o*'! lbs 
m, — 

_ _ 19.7 lbs 

ma g 

_ _ 215.3 lbs 

m3 

_ _ 44.7 lbs 
m. 


Acoustic Power Input 

This term for the external icoustic power input is handled as in Example 

An 

Problems 1 and 2. The applicable term for is 




surface 


= 237,473 Af 


The internal acoustic field was assumed to be reverberant. Therefore, 
the internal acoustic excitation was handled in the same manner as the 
external acoustic excitation. The input term is 

^ _ 2tt Co*A4<Pi>ct N 4 (surface) 

(Oq (Aoj) ni4 
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All of the inodes of this elanent are surface inodes, so 


"'(surface) ’ "»<les-'Hz 

The surface weight density for the cover Is 0.00603 Ib/ln*. Therefore 
the ^ term Is 

IQ 

A (-27 Af) . 17,234 if 


The Internal sound pressure levels were measured during MPTA testing and 
are listed In Table III-2. 

The radiation efficiency values for the DPI box cover are from Appendix II, 
based on a coincidence frequency of 1434 Hz. This coincidence frequency 
for the laminated cover Is given by 



12(l-v2) I 


where the equivalent th1ckn#‘ss for the cover was determined, as In the 
Modal Coupling section of this example problem, to be 

t = 0.189 In. 

The radiation efficiency values are also listed In Table III-2. 
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Table I1I-2 


Frequency 

(Hz)_ 

<^> 

(dB re 2x10*® N/mM 

£ 

50 

122.9 

0.00316 

63 

124.4 

0.00398 

80 

125.6 

0.00501 

100 

126.2 

0.00631 

125 

126.5 

0.00794 

160 

126.5 

0.0126 

200 

126.3 

0.0158 

250 

126.0 

0.0200 

315 

125.5 

0.0251 

400 

124.7 

0.0316 

500 

123.6 

0.03^8 

630 

122.5 

0.0501 

800 

120.8 

0.0631 

lOOC 

119.5 

0.126 

1250 

117.5 

0.501 

1600 

116.5 

3.98 

2000 

115.5 

2.00 
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Response Solution 

The response solutions for each element were determined In each 1/3 octave 
band from 60 to 2000 Hz. A sample prediction for th.. bjrvl.. 4. presented 
below as an example. **** 


(wHi ) El “ Ni^iaEz ~ Si 

•N2<I>12Ei + (usiz Ni<j)i2 + ^34*23 ) E; ~ N24>23E3 ~ 0 
-N34>23E2 (ciVl 3 + N24>:.3 H«»*^3'i ) E 3 - N3<t>3i|E4 * 0 

“Nif^isiiEa + (wn<» H3<I>3‘» ) E 4 - $4 


Substituting for the parameter values, the expressions become 




« r Ro * 

2irH(in6 X 12)"]J237,473) 

|[l0~^ X 8.41 X 10'^* 

[.000794] 

[(2ti X D0)^jf 

2x X 



-['>-32(4’33)£^[T (2, k'sO)'] 

[ 2 . X 50 ][ 0 . 1 ] ♦ [ 6 . 52 (^)][^]» [ 0 . 12 (jf 3 'J^]j 


L 9 ( 2 tt X S 0 )^J I JL ^ ( 2 it X 50 )^J 
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wfiere ^f ~ 4"^ octave bands. 

Dividing by g and solving for the mean squared accelerations; 



= .0135 




.0433 
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The acceleration spectral density levels are 


<a^> 

* •» 

.0135 _ 

g* (Af) 


<i^> 

.130 . 

K> 

-h 

1 

(?^) 

<^> 

.0433 _ 

S* (Af) 

03 

<^> 

.646 , 

g* (Af) 

5b 

03 


.00117 

.0112 

.00375 

.0560 


and the root-inean-squared acceleration in this 1/3 octave band Is 
(g^s)^ -yTnO * .360 

(grms)j * V- 0433 = .208 

(9ms), =VT646 = .804 

The predicted response levels are plotted in Figures III-16 through III-21 
in acceleration spectral density and octave band formats. 

The criterion of 20 modes per analysis band in each element as a requirement 
for SEA prediction accuracy is satisfied tbr the model at 800 Hz, indicating 
predictions with ± 3 dB accu*'acy above this frequency. 
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FREQUENCY »HZ) 


Figure III-16. Predicted Response (g^/Hz) for External Tank Interstage Area 
(Elements 1 and 2 ) 
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Acceleration Spectral Density. g*/Hz 



frequency (hz: 


Figure III-17. Predicted Response (g^/Hz) for External Tank Interstage Area 
(Elements 3 and 4) 
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Figure II 
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Figure III-20. Predicted Response (grms) External Tank Interstage Area - DFI Equipment Panel 
(Element 3) 
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Comparison of Prediction with Measured Test Data 

Subsequent to completion of this analysis, NPTA response data were made 
available for comparison with the predicted responses of elements 1 and 3. 
The external shell (element 1) prediction Is essentially the same as for 
Example Problem Number 1 and has been previously discussed. The data 
comparison for the DPI equipment panel (element 3} Is presented In Fig- 
ure III-22. The comparison shows the prediction to furnish a good 
approximation to the average response at frequencies above 100 Hz. The 
peak In predicted response at 1600 Hz Is due to a coincidence frequency 
effect of the DPI box cover (elanent 4), similar to the external shell 
response peeking discussed In Example Problem Nimiber 1, which was attributed 
to Improper damping definition near the coincidence frequency. 

A supplementary check case was performed to evaluate sensitivity of the 
panel response prediction to variation In the external acoustic levels. 

The external acoustic levels were Increased by 6 dB for the check case 
and resulted In essentially no change In the predicted DPI panel response, 
showing It to be mainly driven by the Internal acoustics. 
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FBfcOUENCY (HZ» 


Figure III-22. Comparison of Measured Test Data (dashed line) with Predicted 
Ttesponse (solid line) for DFI Equipment Panel (Element 3) 
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Example Problem Number 4 
SPACE SHUTTLE RETRIEVABLE SPACECRAFT 


The structure to be analyzed consists of the fairing for a Delta launch 
vehicle and a payload spacecraft. The payload attaches directly to the 
fairing rather than to a lower stage of the Delta launch vehicle. This 
analysis makes use of the averaging abilities of SEA to provide a gross 
estimate of the payload response for evaluation of the attachment configura- 
tion. The analysis utilizes information and data obtained during a number 
of previous modal analyses of the Delta launch vehicle and payload space- 
craft. Since the fairing response to the flight acouslic environment had 
bee’^ previously measured, the system was treated as having a mechanical in- 
put from the fairing to the spacecraft. This treatment assumes no change in 
fairing average response with the spacecraft connected directly to the fairing. 
Since this configuration may be expected to attenuate fairing response some- 
what due to mass loading, the spacecraft response predictions presented herein 
are expected tc be a conservative estimation. 
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Model 


The model for the structure will have two elements, one for the fairing and* 
one for the payload. These elements are illustrated in Figure I I 1-23. 




Figure III-23. SEA Model of Fairing and Payload Spacecraft 
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Response Equations 

The SEA response equations for the two-element system witn external acoustic 
excitation are 


(mHi + N2^12) " Ni4>i2E2 ~ S* 


-N24ii2Ei + (wna + Ni^i2) E 2 = 0 

inhere the subscript 1 denotes external shell v<'1i*es, and 2 denotes DFI box 
and Intercostal values. 

Since the response of element 1 Is known, the response of e1en(>nt 2 can be 
determined using only the second equation: 

E2 * 

I(jn2+I'lj^i2) 


Damping 

A damping value of 1-1/2% (n=0.03) was used for both of the model elements. 
This damping value had previously been selected for use with modal analyses 
of the Delta fairing and payloads. 


Modal Density 

The modal density of the fairing was calculated by assuming the isogrid 
structure to be composed of beam and plate elements. Modal densities 
were then calculated for each of the beam and plate sub-elements, and these 
sub-element values were summed to give the total modal density of the 
fairing. 

The modal density of a single plate sub-elanent is: 


np(<.)=. 


4ith 


12p(l-x>^) 
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and for a beam: 


n. (tit) = L 

Totalling the modal densities: 


n(w) = n. + n = 1 

V-v 1 

\ ^ 

^ 

y ^ 



/ E 

Z h 



4ff V12 p(1-'»‘^ 



The summations account for all the modal properties of the stiffeners and 
skin. The values for these siii. used for the fairing were: 

Tm = ^33. VI. in S~ = 2,181.000 in 

Z-h 

These modal densities were corrected for curvatu. e effects based on a ring 
frequency of 650 Hz. erection factors were obtained from Reference 4. The 
resulting number of modes for each frequency band of analysis is listed in 
Table III-3. The plate medal density equation presented is actually for free- 
free boundary conditions, hul yields a valid approximation for all boundary 
conditions at high frequencies. However, this technique is expected to 
overestimate the modal density in the lower frequencies of analysis. 

The modal density of a "typical" peyload is estimated to be approximately 
the 3 amc as for the fairing. This estimate is based n comparisons of 
the results of modal analyses for the fairing end for 2000-pound and 4000- 
pound Delta payloads. The results of the ana’^yses, presented in Figures 
III-24 and III-25, yielu essentially ide.utical estimates for the modal densities 
of the two elements. Although simple, beam-type models were used for 
* ese enalys^., the relat've equality of tne modal densities may be expected 
io remain when the higher fi-eouency modes are ccnside.vd. 



Table III-3 


f 

N 

PSD 

(sVhz) 

250 

633 

3.0 

315 

990 

4.5 

40G 

1500 

15.0 

500 

2320 

18.0 

630 

3960 

10.0 

800 

4780 

3.5 

1000 

5280 

0.8 

1250 

6230 

0.5 

16U0 

7410 

0.45 

2000 

9120 

0.45 
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Modal Coupling 


In order to estimate the coupling parameter, the structural joint between 
the Saturn S-II aft skirt and thrust structure assemblies was selected as 
typical of the fairing/payload joint. The S-II aft skirt/thrust structure wa . 
modeled as an SEA system of two elements with only tne aft skirt externally 
excited. This model permitted solving for the coupling parameter of the 
system using the explicit response solution 



2 

<3j > 


A 



Nj* 

wht 


where "T" subscripts indicate the thrust structure and "A" subscripts indicate 
the aft skirt. 


Assuming thrust structure damping to be a maximum of 1 percent and that the 
coupling decreases with frequency at the same rate which MDAC experienced 
with the UpSTAGE evaluation, the S-II aft skirt/thrust structure SEA coupling 
parameter was estimated to have the values shown in Figure III-26. These 
values were used for the coupling parameter of the fairing/payload model. 

Element Energy 

The element energy was handled as in Example Problem 1. 


1348 lbs 

mi 

_ 390Q lbs 
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Figure III-26. Field-Joint Coupling (#(,) for Delta/Payload Analysis 
(Based on S-II Data) 


The measured response of elonent 1 is shown in Figure IIi-27 in an accelera- 
tion spectral density format with units of 


PSD 


9"^ (Af) 


Therefore 


<ii;> = g2(Af)*PS0 

C3i! be evaluated from Figure 1 1 1-27 by using PSD value for the center 
frequency as representative over the bandwidth, Af. The PSD values used for 
analysis are listed in Table III-3. ' 
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Figure III-27. Fairing Skin Response (Radia’i) 


Re sponse Solution 

The response predictions for each element were determined in each one- third 
octuve band from 250 to 2000 Hz. A sample prediction for the 250 Hz band is 
presented below as an example: 


(una + N, (^17) 


Substituting for the parameter values, the expression becomes 

(6.33)(1.75) 


3900 <&l> 


g ( 27 tX 250 )? |[ 2 T 7 X 250 ]!; 0 . 033 +[( 633 )( 1 . 75 )]j 


/250_ 

1348 g-\4. 33/(3. 0) 


9 (2^X250)^ 


where Af = . 4~33 one- third octave bands. 
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The acceleration spectrcl density level is 


-a|> 


g^(Af) 


0,43 


and the rool-mean-squared acceleration in ^his one-third octave band is 


/<a2> 




The predicted response levels are plotted in Figures III-28 and III-29 
in acceleration spectral density and in one- third octave band formats. 

As previously noted, these levels can be expected to represent a conservative 
overestimate of the actual response due to an expected attenuation of the 
fair response in this configuration. 
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Acceleration Spectral Density, g^/Hz 



FREQUENCY (HZ) 


Figure III-38. Predicted Response (g^/Hz) for Payload Element 
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Appendix IV 


SEA RESPONSE SOLUTION PROGRAMS FOR 
HEWLETT-PACKARD AND TEXAS INSTRUMENTS CALCULATORS 


SOLUTION FOR MODEL RESPONSES 

One of the final steps of the SEA applications procedure will be the mathe- 
matical solution for the system response. Because of the repetitive nature 
of the response calculations for the various frequency bands of an analysis 
(often the 14 one-third octave bands fran 100 to 2000 Hz will be Included 
In an analysis), response calculations can be expediently accomplished through 
the use of preprogramned solutions. Such programs also eliminate the errors 
which occur with monotonous, repetitive, hand calculations. The capability 
of many programmable hand calculators Is suitable for solution of the smaller 
conmon model sizes. Programs for both Hewlett-Packard and Texas Instruments 
programmable calculators are provided, thus encompassing the more popular 
equipment In current use, as well as providing examples In both reverse 
polish and algebraic notation for modification to additional calculator 
systems. 

HEWLETT-PACKARD (HP-67, HP-97) SEA RESPONSE PROGRAM - TWO ELEMENTS 
Description 

This program provides the SEA response solution for a system of two elements. 
The response equations are 

(wrii + N 2 (|>i 2 ) Ex - Ni4)i 2E2 = Sx 
-N2<|)X2Ex + (wn2 + Nx<|)x 2 ) E 2 “ S 2 

The energy Input terms may be input directly to the program or can be cal- 
culated if a reverberant acr .stic fielu provides the structural excitation. 

The required Inputs to the program are listed below: 
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f - center frequency of analysis band (Hz) 

Af - bandwidth (Hz) 

ni - element 1 damping factor (dimensionless) 
j)z - element 2 damping factor (dimensionless) 

<)>i 2 - coupling value for elements 1-2 sec ') 

ni - element 1 modal density (modes/Hz) 
n 2 - element 2 modal density (modes/Hz) 

Wj - weight of element 1 (lbs) 

Vl 2 - weight of element 2 (lbs) 

Si - element 1 energy input ) 

$2 - element 2 energy input 


Additional requirements if the are to be calculated for a reverberant 
acoustic field; 

Cq - speed of sound in surrounding medium (’n/sec) 

(^) ‘ ^ surface mass parameter 

(^) " ^ surface mass parameter ( ]b-s ' ec2- H z ) 

(SPL)j - element 1 sound pressure level (dB) 

(SPL)^ - element 2 sound pressure level (dB) 

02 - element 1 radiation efficiency (dimensionless) 

02 - element 2 radiation efficiency (dimensionless) 

Output from the program is the root -mean- squared acceleration of the elements 
over the input frequency interval; 


<^> 




and 




g" 
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USER INSTRUCTIONS 


Step 


Enter 

Press Display 

0 


— 

RTN 

— 

1 

Input number of elements (N=2) 

2.0 

R/S 

2.0 

2 

Input center frequency 

f 

R/S 

2irf 

3 

Input bandwidth 

Af 

R/S 

Af 

4 

Input damping of elanent 1 

Hi 

R/S 

Hi 

5 

Input damping of elanent 2 

r \2 

R/S 

Hz 

6 

Input element coupling 


R/S 

•^12 

7 

Input modal density of element 1 

ni 

R/S 

Ni 

8 

Input modal density of elanent 2 

Hi 

R/S 

Nz 

9 

Input weight of element 1 

Wi 

R/S 

Wi 

10 

Input weight if elanent 2 

Wz 

R/S 

Wz 

11s 

Input "0" If element energy input 
terras (S^) are to be Input to 
program; Input "1" If reverberant 
acoustic Input Is to be calculated 
by program and see below. 

0 

R/S 

0 

12a 

Input Si 

Si 

R/S 

Si 

13a 

Input Sz, obtain solution 

Si 

R/S , 

V 

pv 

' 9 ^ 

14a 

Obtain solution for e'-t.t - t 2 

/<if> 

V 

R/S J 

j<g> 

9 ^ 

15a 

Check for end of program 


R/S 

0 

•lb 

Alternate solution with calcula- 
tions for reve»-berant acoustics 

1 

R/S 

1 

12b 

Input speed of sound (in/sec) 

Co 

R/S 

Co 

13b 

Input (“) fo • element 1 

(^). 

R/S 


14b 

Input for element .? 


R/S 

(An) 

\ m /2 

15b 

Input SPL for elanent 1 

(SPL)i 

R/S 

(S. )i 

16b 

Input SPL for el .ament 2 

(SPL) 2 

R/S 

v3PL)2 

17b 

Input radiation efficiency for 
element 1 

Ol 

R/S 

Oi 

lUb 

Input radiation efficiency for 
element 2, obtain solution 

02 

R/S J 

9^ 


90 





Enter 

Press 

Display 

19b 

Obtain solution for element 2 

J? 

R/S 

p>" 

20b 

Check for end of program 

— 

R/S 

0 


Data Registe* 




A 

ni 

7 

a 12 

B 

H2 

8 

a?- 

C 

Wi 

9 

322 

D 

Wa 

PO 

— 

E 

ni 

PI 

Cq 

I 

(27Tf)V386 

P2 

Oi 

0 

na 

P3 

Oz 

1 

<l>12 

; 

(SPL)i 

2 

Si 

1 P5 

(SK)2 

3 

S2 

P6 

(^) 

4 

N ] 



5 

2TTf 

P7 

\m Iz 

6 

3ll 1 

' P8 




! P9 

... 


where the ai'j represent program-generated matrix elements. 
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rro^ran 


Progran 
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- 
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.C »l - 


' ■ ~ — 

. 

115 . 

<at 


SIS'* 


45? 

* 

U6 

s 

Cl 5 

a/s 


363 

?sS 

ij7 

STS5 




361 

iiu,5 

lie 

aci£ 

tfC*? 

A 


36i 

k* 

115 

RCt,i 

?tf6 

2 • 


36S 

4 

iZ3 

X 

9tf“ 

Y 


fiS« 

PSs 

Ui 

ROS 

m 

5755 


««5 

STM 

l££ 

X 

ciJ 

S' 5 


m 

3CL6 

u3 

ROJ 


S7Ji 


367 

A 

124 

* 

• * « 

vi 


368 . 

S8Z2 

i2f 

CHS 

ii; 

iT&t 


dS| ? 

JC . ' ” 

U6 

3TS7 

41^ 



ST2 

3CL4 

127 

803 

ai4 

SICS 

V - 

«7l 

1 

:2S 

RCLi 

m 

i'i 



3 

’ 125 

A 

iio 

JTii 


673 

i 


xcu 

HT 



37* 

134 

iJi 

xi 

ii6 

«cu 


drs 

X 

152 

SOI 

dii 

X 


376 

FSS 

155 

* 


3TC£ 


arr 

STC2 

234 


«e: 

i'^5 


0?S 

KS 

155 

S766 

ili 

RCLi 


373 

scl: 

IIo 

ROS 

«e? 

A>* 


333 

3U7 

IT 

RCLB 

€£4 

ST53 


36: 

X 

iSs 

X 

«15 

a^s 


3SZ 

3CL5 

139 

RCLc 


sm 


335 

X 

243 

RCLi 

«£• 

ft/£ 


384 

sa.5 

141 

w 

Jh 

ezi 

ST6L 


365 

1 

142 

T 

6ZS 

«/3 - - 


386 

—e - . . - - 

145 

RCL5 

333 

x*«-? 


367 

♦ 

1*4 

X 

vSl 

S7S1 


363 

16“ 

145 

RCLI 

£•• 

v««» 

rC£ 


383 

* 

146 

* 


«'S 


353 

ftS 

147 

STOS 

334 

sm 


891 

sro3 

1*8 

1/X 

335 

8^S 


356 

S195 

149 

SCwS 

ea 

ST3£ 


857 

4LoLl 

158 

X 

337 

R'S 


354 

R^S 

151 

RCm7 

333 

ST67 


355 

STM 

*e; 

X 

333 

i'S 


656 

RS 

in 

CHI 

643 

STC4 


357 

STC3 

154 

ROl 

3^1 

X^S 


398 

kM2 

155 

t 

3*i 

ST35 


695 

2£L5 

156 

RClS 

343 

2/0 


136 

hi 

1ST 

RCL9 

344 

S'ii 


I3i 


158 

4 

tf-5 

r/£ 


i62 

8 

l5i 

RCl7 

34i 

STvS 


135 

6 

16*^ 

X 

34T 

f; 


134 

4 

161 

Citi 

343 

acLi 


i65 

sTc: 

162 

X06 

343 

X* 


136 

PCt.5 

:65 

T 

353 

y 


787 

nCLh 

1£4 

* 

c5. 

8 


168 

X 

:£5 

S?w 

C52 

• 


:35 

Ki^m8 

:$0 

irC 

353 

4 


113 

RCLI 

167 

7X 

354 

1 

• 


• * • 

*4 * 

X 

168 

«TA 

355 

££a 


4*«» 

t 

s»9 

R'$ 

S5i 

i 


IIJ 

xu: 

1/9 

.r* 

45? 

a 

c 


X 

92 

1 ♦! 
•f ’ 
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Progros 


Step 


172 

- 

» 

jTj 



RCl3 

i Tr 

♦ 

n 


• • 

ifc 

« A 

17? 

pf>rA 

Jic 

R'S 

is: 

V 

161 

PTs 

ill 

hi 


TEXAS INSTRUMENTS (TI PROGRAmABLE 59) SEA RESPONSE PROGRAM - TNO OR THREE 
aEMENTS 

Description 

This program provides the SEA response solution for a system of two or three 
elonents. The response equations are: 


(wn» N2 $‘i 2 ■♦’“Na^ia) Ei - N14112E2 - Ni^^iaEa * Si 


•N2^i2Ei ■*’(u>n2 Ni(J»i 2 ^ Na<^2a) E2 ~ N2<^2aE3 * S2 
“Na(j>iaEi ~ Nai|>23E2 +(wna + Ni4>.3 + Nat^aa) E3 * Sa 


The energy Inpu^ terms. S{, may be input directly to the program or can be 
calculated If a reverberant acoustic field provides tte structural excitation. 
The required Inputs to the program are listed below: 

f - center frequency of analysis band (Hz) 

Af - bandwidth (Hz) 

m - element 1 damping factor (dimensionless) 

H 2 - element 2 dz.TipIng factor (dimensionless) 

Ha *• element 3 damping factor (dimensionless) 
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♦l* " 

♦l3 ' • - 
♦*3 “ 

ni - 
nz - 

n$ - 

Ml - 

M* - 

Ml - 
Si - 

$2 - 

$9 - 


coupling value for elements 1-2 

j- 

coupling value for el^nts 1-3 ( poa^gi^ ) 

coupling value for elements 2-3 ( 55 ^: 15 ^) 

element 1 nndal density (modes/Hz) 
element 2 modal density (modes/Hz) 
elonent 3 modal (tensity (modes/Hz) 
weight of element 1 (lbs) 
weight of element 2 (lbs) 
weight of element 3 (lbs) 
eien»nt 1 energy Input ( 

element 2 energy input 

element 3 energy Input 


Additional requirements If the are to be calculated for a reverberant 
acoustic field: 


Co 

m 

(^)3 

(SPL)i 

(SPL)i 

(SPL)s 

Oi 

Oz 


03 


speed of sound in surrounding medium (In/sec) 

/in*-roodes \ 

element 1 surface mass parameter lib-jec^-Hz' 
elanent 2 surface n»ss parameter 

/in*-modes \ 

element 3 surface mass parameter \ iblsec^Hz / 

element 1 sound pressure level (dB) 
element 2 sound pressure level (dB) 
element 3 sound pressure level (dB) 
element 1 radiation efficiency (dimensionless) 
elanent 2 radiation efficiency (dimensionless) 
element 3 radiation efficiency (dimensionless) 
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USER INSTRUCTIORS 


Procedure « W « 2 (see following procedure for N»3) 

Step Enter Press Display 


Oa Partition oenory < 

Ob Insert progrMi 

Oc Initialize program 

1 Input nunber of elements* N 2 

2 Input center frequency f 

3 Input ban(to1dth Af 

4 Input dancing of element 1 m 

5 Input dai^tlng of elo^t 2 , na 

6 Input element coupling 4is 

' T' Input modal dafisity of element 1 ni 

8 Input modal density of element 2 nz 

9 Input weight of element 1 Mi 

10 Input weight of elanent 2 Ha 

11a Input "0" If element energy Input 0 

terms are to be Input to program; 

Input "1“ If reverberant acoustic 
Input Is to be calculated by- 
program* and see below 

12a Input Sj 

13a Input Sa, obtain solution Sa 


• 14a ' Obtain solution for e1en»nt 2 
15a Check for end of program 


Op 17 639.39 


RST 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 

R/S 


Prevlims 

t-register 

2wf 

Af 

ni 

2 

2 

ni 

2 

Ml 

386 

0 



11b 

Alternate solution with calcula- 
tions for reverberant acoustics 

1 

R/S 

0 

12b 

Input speed of sound (1n/sec) 

Co 

R/S 

Co 

13b 

Input (^) for element 1 

(^) 
N m /i 

R/S 

(t). 

14b 

Input for element 2 

(-) 

\ m fz 

R/S 

2 

15b 

Input SPL for element 1 

(SPL), 

R/S 

(SPL), 

16b 

Input SPL for elonent ? 

(SPL). 

ft/s 

2 
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StM 


Enter 

Eim 

§iSfi!3£ 

17b 

Ii^Nit radiation efflclenqr for 
elanent 1 


R/S 

di 

18b 

InfMt radiation efflclenqr for 
elanent 2» obtain solution 

02 

R/S 

i^> 

1 9* 

19b 

(H>ta1n solution for elanent 2 

— 

R/S 

l¥ 

20b 

Check for end of prograa 


R/S 

0 


ProcfiAire ~ H«3 


Oa 

Partition library 

4 

Op 17 

639.39 

Ob 

iisa t program - 

— 

— 

— 

Oc 

Initialize progr^ 


RST 


1 

ItqMt nund)er of elemaits 

3 

R/S 

Previous 

t-register 

2 

Input center frequency 

f 

R/S 

2irf 

3 

liqjut bandwidth 

Af 

R/S 

Af 

4 

Irqmt dashing of elanent 1 

ni 

R/S 

ni 

5 

Input damping of element 2 

H2 

R/S 

2 

6 

Input damping of element 3 

na 

R/S 

ns 

7 

Irq)ut coupling for elanents 1-2 

^12 

R/S 

2 

8 

Input coupling for elements 1-3 

<^13 

R/S 

4i3 

9 

Input coupling for elanents 2-3 

^23 

R/S 

423 

10 

Input modal density of element 1 


R/S 

"l 

11 

Input modal density of elanent 2 

Oj 

R/S 

2 

12 

Input modal density of element 3 

Os 

R/S 

"3 

13 

Input weight of element 1 

w, 

R/S 

Ml 

14 

Input weight of element 2 

^2 

R/S 

2 

15 

Input weight of element 3 

W, 

R/S 

386 

16a 

Input "0" If element energy Input 
terms are to be Input to program; 
input "1" If reverberant acoustic 
Input Is to be calculated by 
program and see below 

0 

R/S 

0 
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Press 

Display 

17a 

Input S| 


R/S 

Si 

18a 

Input $2 

, S 2 

R/S 

2 

19a 

Input S 3 » obtain solution 

S, 

R/S 

/ 

V 9 * 

20 a 

Obtain solution for elo^t 2 

/ 

R/S 

IW 

21a 

Obtain solution for element 3 

— 

R/S 


22 a 

Check for end of program 

— 

R/S 

0 

16b 

Alternate solution with calcula- 
tions for reverberant acoustics 

1 

R/S 

0 

17b 

Input sf^ed of sound (In/sec) 

Co 

R/S 

Co 

18b 

Input for element 1 

(n 

R/S 
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\ ro /j 
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Input for element 2 
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R/S 
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R/S 
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R/S 
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22 b 

Input SPL for element 2 

(SPL ) 2 

R/S 

2 

23b 

Input SPL for elanent 3 

(SPL ) 3 

R/S 

(SPL) 3 

24b 

Input radiation efficiency for 
element 1 

Oi 

R/S 

Ol 

25b 

Input radiation efficiency for 
element 2 

02 

R/S 

2 

26b 

Input radiation efficiency for 
elanent 3, obtain solution 

Oi 

R/S 

1 <«?> 
i 9^ 

27b 

Obtain solution for element 2 

— 

R/S 

V 

28b 

Obtain solution for element 3 

— 

R/S 

<^> 

29b 

Check for end of program 


R/S 

0 
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12 
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ni 
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19 

ni 




where the represent program-generated matrix elanents. 
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